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MULTISTAGE COMPRESSOR SIMULATION APPLIED TO THE 
PREDICTION OF AXIAL FLOW INSTABILITIES 
by Ross G. Willoh and Kurt Seldner 
Lewis Research Center 

SUMMARY 

A simulation model for the one-dimensional steady state and dynamic performance 
of an axial flow compressor is presented. The modeling technique is applied to the eight 
stage J 85-13 compressor. It is shown that the simulation stability boundary coincides 
with the compressor stall line. Analytical and experimental results are presented and 
compared. 


INTRODUCTION 

The transient performance of aircraft propulsion systems, using axial compressors, 
is strongly influenced by the compressor stall characteristic. Compressor stall is ac- 
companied by reduced airflow and engine thrust and increased turbine inlet temperature. 
In supersonic applications the reduced airflow can result in the inlet unstarting. The 
ability to adequately predict stall and to evaluate the effects of operational parameters on 
engine stall characteristics is thus of considerable significance in propulsion system 
design. 

As indicated in reference 1, there is a minimum flow point for each compressor 
operating speed, where the compressor will stall. The locus of minimum flow points 
defines a compressor stall line. Jet engine fuel controls schedule propulsion system pa- 
rameters to maintain engine operation along an operating line which is below the stall 
line. The separation between the stall and steady state operating lines represents engine 
stall margin. 

During acceleration from one steady state condition to another, engine operation fol- 
lows a path from the steady state operating line toward the stall line and back. If the ac- 
celeration time is to be minimized, the stall margin during the transient must be small. 
Knowledge of the compressor stall line is thus of great importance in the design of pro- 



pulsion system controls. 

Much experimental and analytical work has been done on compressor stall. Refer- 
ence 1 contains a compilation of some of this material. A discussion of the influences of 
compressor stall on fuel control design is presented in reference 2. This report presents 
the results obtained from the application of an analog simulation of a compressor and its 
gas dynamics to the prediction of compressor stall. 

The simulation techniques used in this report are similar to those of reference 3 in 
that both were based on the use of stage stacking and lumped volume gas dynamics. The 
techniques of this report, however, use both stage pressure and temperature rise char- 
acteristics and avoid the problems of differentiation through the application of a momen- 
tum balance to determine stage flow. Similar techniques are reported in reference 4. 

The simulation was originally prepared to support an experimental investigation of 
airbreathing propulsion system dynamics and control. The experimental pregram utilized 
the J 85-13 engine, and, hence, the analytical model is based on this engine. The tech- 
niques used, however, can be applied to other airbreathing engines. It was found that the 
simulated J 85-13 compressor exhibited stability problems in the vicinity of the actual 
compressor stall line. Specifically, the stability boundary of the simulated compressor 
matches the actual compressor stall line quite closely. 


MATHEMATICAL ANALYSIS 

The analysis of multistage axial flow compressors presented in this report is based 
on the representation of individual stages by steady state stage characteristics and by the 
gas dynamics of the associated stage volume. An idealized stage, of a multistage com- 
pressor, is shown schematically in figure 1. The symbols used in this figure and through- 


(n-l) stage — *--* n stage •— «— (n+1) stage 


I 1 1 1 



Figure 1. - Schematic of n** 1 compressor stage. 
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out the text are defined in appendix A. 

The analysis assumes that each stage, consisting of a rotor, a stator, and an included 
volume, can be modeled by a compressor element coupled to the gas dynamics of a lumped 
volume element. The compressor element is assumed to instantaneously satisfy the stage 
steady-state characteristics in response to the conditions imposed by the adjacent vol- 
umes. The details of the mathematical representation of the steady-state stage perform- 
ance are presented in appendix B, and the equations for the lumped volume gas dynamics 
are developed in appendix C. 


Steady-State Characteristics 

The prediction of the off design steady- state performance of a multistage axial flow 
compressor is a complex problem. The inclusion of gas dynamics, required to examine 
system stability, further complicates the problem. Reference 1 includes discussions of 
the blade element method, stage-stacking techniques, and the use of performance data 
from existing compressors for the prediction of off design compressor performance. Of 
the three techniques, only the stage- stacking method combines easily with lumped volume 
dynamics to form an overall dynamic compressor model. The analysis of compressor dy- 
namics presented in this report is thus based on the stage stacking method for the steady- 
state portion of the model. 

The steady-state performance of the individual stages of a multistage compressor 

P T 

can be represented by pressure and temperature rise coefficients i// and , which 

V C 

are functions of the flow coefficient <f> . These coefficients are mathematically defined 

iu C 

for the n- stage in appendix B as 
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It is shown in appendix B that the axial velocity v z can be related to the stage inlet static 
pressure, static temperature, and flow by 

W T „ JR 

v zc n = ' ■- y) 1 (6) 

n A P 

C, n sv, n-1 

nr p 

The coefficients i// and i// of equations (1) and (2), however, are defined in terms of 
total pressures and temperatures. An equation relating v to stage inlet total pressure 
and temperature is developed in appendix B as 


W \]e ~T v 
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where 6 and 0 are total pressure and temperature ratios defined by 
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The angle j3, defined as the angle between the air velocity and the axial direction, is 
primarily a function of compressor design and speed of rotation. 
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Dynamic Characteristics 


The representation of compressor gas dynamics is based on the assumption that 
overall compressor dynamics can be simulated by individual lumped volume dynamics for 
each stage. The individual stage dynamics are assumed to occur in an equivalent stage 
volume which is located just downstream of the stage's active element. The volume dy- 
namics are modeled with conventional lumped volume techniques. 

A schematic drawing of a typical stage with its associated volume is shown in fig- 

tVi 

ure 1. Through the application of continuity, momentum, and energy balances to the n 
stage, it is shown in appendix C that 
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( 10 ) 


The equation of state for the gas in the volume can be written as 

\5/2 


/ 9 \°/ z 

P tv, n = R ( 1+0 - 2 <n) (P S v,n T tv,n> 


( 11 ) 


Equations (8) to (11) form a complete set of lumped volume dynamic relations for the n^ 1 
stage volume. 


Computer Solution 

Equations (1) to (5) and (7) to (11), when combined with the appropriate stage maps, 
form a set for the n th compressor stage. This set provides a mathematical representa- 
tion for both the dynamic and static characteristics of the n stage of a multistage axial 
compressor. The computer mechanization of the equations can be accomplished with 
either analog or digital techniques. The analysis of this report, however, is based pri- 
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marily on the results of an analog computer simulation. The discussion of the computa- 
tional techniques will thus be directed primarily toward an analog representation. 

The steady-state compressor representation utilizes plots of the normalized pressure 
coefficient if/ and the normalized temperature coefficient \fs plotted against the nor- 
malized flow coefficient $. Figure 2 contains plots of the pressure and temperature co- 
efficients against the flow coefficient for the second stage of the compressor studied. From 
these curves it can be seen that, for a given stage pressure rise the flow coefficient $ is, 
in general, double valued. It is thus necessary to enter the stage maps with the flow 
function. 

Equations (3) and (5) can be used to compute the flow coefficient $ from the stage 
inlet weight flow, static pressure and temperature, and the compressor rotational speed. 
Equation (3) can, however, be rewritten in the form 



( 12 ) 


The term v z n /V©^ of equation (12) can then be computed from equation (7). This rep- 
resentation allows the use of total pressure and temperature in computing the flow coef- 




Figure 2. - Pressure and temperature maps for 
second stage of compressor. 
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Figure 3. - Air velocity computation. 

ficient <&. 

The direct solution of equation (7) for each compressor stage would severely compli- 
cate an analog simulation. The axial velocity term v z /V®^ computed from equation (7) 
for various corrected weight flows and rotor air inlet angles is plotted in figure 3. The 
range of corrected weight flow and rotor air inlet angles included in figure 3 covers the 
compressor normal operating region. The deviation from a mean value of air inlet angle 
for stages 2 to 8 of the compressor are plotted as functions of compressor corrected 
speed N in figure 4. The largest rotor air inlet angle encountered in operation is 
less than 26°. It can be seen from figure 4 that, with the exception of the second stage, 
none of the angles vary more than ±2° from the mean value as the corrected speed varies 
from 60 to 100 percent. 

It can be seen from figure 3 that the axial velocity to corrected flow characteristic 
is not extremely sensitive to variations in the air inlet angle. The variation of ±2° or 
less in the air angle for stages 3 to 8 was thus neglected and the stages were assumed to 
have constant rotor air inlet angles. The second stage air inlet angle varies o|° over the 
60 to 100 percent speed range. As shown in figure 3, however, the axial velocity to cor- 
rected flow characteristic is much less sensitive to changes in the air inlet angle at low 
mean values of the angle than it is for high mean values. The second stage, which has 
the lowest mean value of rotor air inlet angle, was also assumed to have a constant rotor 
air inlet angle. 
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Corrected speed, percent 

Figure 4. - Deviation of stage rotor air inlet 
angle from mean. 


For the purposes of the computer representation, the axial velocity to weight flow 
characteristic for each of the stages 2 to 8 was based on a constant mean rotor air inlet 
angle. The actual simulation used a quadratic curve fitted to the appropriate character- 
istic. 

The J 85-13 compressor contains variable inlet guide vanes ahead of the first com- 
pressor stage and variable interstage bleed in the third, fourth, and fifth stage volumes. 
The two variable geometry features are scheduled, between their extreme limits, as a 
function of engine corrected speed N /V®* The variation occurs in the corrected speed 
range from 80 to 100 percent. 

The effects of the variable inlet geometry are included by varying the first stage 
pressure and temperature coefficients \p c j and if/* ^ as functions of corrected engine 
speed. Computationally, the variation was achieved by generating the pressure and tem- 
perature coefficients as functions of the flow function for the two extremes and an inter- 
mediate case and interpolating linearly between the three values for each function. 
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Although the simulation of the individual stage gas dynamics is essentially based on 
equations (8) to (11), modifications were included for three of the stages to account for 
interstage bleed. Further, the distinction between total and static quantities was dropped 
in the computer representation of the stage dynamics. This simplification, which is dis- 
cussed in appendix C, is essentially equivalent to eliminating the Mach number terms 
from equations (9) and (11). With this simplification these equations become 


dt 


(W 


c, n 


A„g 

)=— (P, 


p ) 

tc, n tv, n ; 
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(13) 
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(14) 


The interstage bleed is computed for the third, fourth, and fifth stages by assuming 
the bleed ports are choked and that the bleed flow can thus be computed from 


^b, n k b A b, n pjr — 
V 1 tv, n 


(15) 


The bleed port area n is determined from engine operating schedules. The effects 
of the bleed flow are included in the stage representation by modifying the continuity bal- 
ance of equation (8) and the energy balance of equation (10). For the stages with bleed, 
these equations become 


j(Psv,n)=^<W c ,„-W c , n+ 1 -W b;n ) 


(16) 


^ ^sv, n^tv ? ~ y < T tc, rftc, n " T tv, n^c, n+1 ~ T tv, n^b, 


(17) 


n 


A block diagram illustrating the computational sequence used for a typical stage is 
presented in figure 5. The stage inlet conditions are combined to determine the flow 
function. The flow function is used to determine the stage pressure and temperature co- 
efficients, and, hence, the stage discharge total pressure and temperature. The dis- 
charge conditions are then used with the lumped volume gas dynamics to compute the 
stage volume gas properties, and, hence, the inlet conditions for the succeeding stage. 

The overall J 85-13 compressor representation was formed by interconnecting indi- 
vidual stage representations for each of the eight stages. The first stage maps were 
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Figure 5. - Block diagram of simulation for n* h compressor stage. 


modified, as discussed, to account for the variable inlet guide vanes, and bleed flow ef- 
fects were included in the third, fourth, and fifth stages. 


RESULTS AND DISCUSSION 

The stage stacking technique used for the representation of the compressor steady- 
state characteristic requires the simulation of individual stage characteristics for each 
stage. To simplify the initial simulation, however, it was assumed that the dynamics of 
several stages could be combined and represented as a single lumped volume element. 

The stage volumes with interstage bleed ports were individually simulated to maintain the 
bleed flow dynamics. The dynamics of stages 1 and 2 and the dynamics of stages 6 to 8 
were thus combined to form two lumped volume elements. 

The compressor simulation was terminated with an equivalent combustor volume and 
a choked exit station, and the compressor variable geometry was maintained on normal 
operating schedules. The simulated compressor was then mapped by varying the choked 
station area at various rotational speeds. With increasing compressor terminating im- 
pedance, caused by reductions in the choked station area, the compressor operating point 
moves along a constant speed line. Figure 6 presents, in schematic form, a typical con- 
stant speed operating point locus. 

It was found that, as the terminating impedance was increased, the simulation would, 
at some operating point along the constant speed line, begin to exhibit a limit cycle. Fur- 
ther increases in the impedance resulted in increased limit cycle amplitude. Typical 
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Figure 6. - Schematic of simulation constant speed 
line and stability boundary. 


limit cycle amplitudes are also indicated in figure 6, The limit cycle exhibited hysteresis 
in that, once the oscillation was established, it was necessary to reduce the terminal im- 
pedance beyond the point at which the oscillation first occurred before stable system 
operation could be reestablished. 

A number of constant speed lines were mapped. Noting the operating condition on 
each speed line where the self- sustained oscillations first occurred allowed the determin- 
ation of a stability boundary for the compressor simulation. This stability boundary is 
plotted in figure 7. For reference, the normal stall line for the compressor is included 
in the figure. The stability boundary resembles, but does not match, the compressor 
stall line. 

Noting that the overall compressor simulation is quite nonlinear, it was concluded 
that the volume lumping used in the initial simulation was not justified. The simulation 
was thus revised to include a lumped volume dynamic element for each compressor stage. 
The stability boundary for the revised simulation was obtained and is also plotted in fig- 
ure 7. It can be seen that the stability boundary of the compressor simulation falls quite 
close to the actual compressor stall line when individual stage dynamics are included in 
the computer representation. 

Since the actual compressor stall line, used in the initial portion of the study, was 
based on minimal data, the technique was applied to operating conditions for which test 
data was available. Compressor performance and stall line data obtained from tests per- 
formed on a J 85-13 engine is presented in reference 5. The experimental results pre- 
sented in this reference are not representative of standard engine operation. The test re- 
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12 14 16 18 20 

Corrected airflow, kg/sec 

Figure 7. - Compressor stall lines determined from a 
grouped stage volume simulation and from an indi- 
vidual volume simulation. 


suits were obtained with the first stage turbine nozzle area reduced 26 percent, and with 
the engine variable geometry operated off the normal design schedules. 

The tests were performed by maintaining constant engine speed through manual ad- 
justment of the power lever, while the exhaust nozzle area was reduced in small incre- 
ments. The equilibrium point just prior to compressor stall was defined as the stall 
point. The procedure was repeated at various engine speeds to obtain the compressor 
stall line. 

The stability boundary for the analog compressor simulation was obtained for the 
same operating conditions as the test data. The simulated variable geometry was sched- 
uled to match the test conditions, and the terminating impedance was varied while the 
compressor rotational speed was held constant. Constant speed lines were plotted and 
the stability boundary was noted. The test data from reference 5 and the results obtained 
from the simulation with individual stage volume representations are both plotted in fig- 
ure 8, The good correlation between the experimental and analytical results can be seen 
in this figure. 

Some insight into the overall operation of the compressor in the vicinity of the stabil- 
ity boundary can be gained from examining the normal compressor operating point and 
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Corrected airflow, kg/sec 


Figure 8. - Compressor map for off schedule inlet guide vanes and interstage 
bleeds. 


Stage 



•50 .55 .60 .65 .70 .75 .80 .85 

Normalized flow coefficient, $/$ 

Figure 9. - Stage pressure maps including operating points for design speed. 



an operating point just prior to the limit cycle. Figure 9 contains the pressure coefficient 
curves for each stage and the two operating points at 100 percent corrected speed. Fig- 
ure 10 contains the same data obtained at 80 percent corrected speed. The compressor 
variable geometry is fully closed at 100 percent corrected speed and fully open at 80 per- 
cent. The effects of the variable geometry are reflected in the different first stage char- 
acteristics and in the reduced flow through the later stages at the 80 percent operating 
condition. 



Figure 10. - Stage pressure maps including operating points for 80 percent of design speed. 

From the data plotted for the 100 percent corrected speed case in figure 9, it can be 
seen that the normal operating condition falls on the negative slope portion of the pressure 
coefficient curves for all of the stages. As the stability boundary is approached, all of 
the stages move toward operation at a lower negative slope. Second stage operation falls 
on the horizontal portion of the stage characteristic and the third stage is operating in a 
region with positive slope. For operation at 80 percent corrected speed shown on fig- 
ure 10, however, normal operation for the second and third stages is always in the region 
of positive slope. As the stability boundary is approached, operation of these two stages 
shifts toward increasing positive slope. 

The stage characteristics plotted in these figures are those included in the simulation. 
It should be noted that they include no sharp discontinuity to represent stage stall. The 
limit cycle encountered in the simulation, that matches the compressor stall line, is thus 
assumed to be a system stability phenomenon associated with the system gains and phase 
lags. Stall prediction by this method cannot be expected to duplicate stall caused by ex- 
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treme degradation of individual stage performance at a specific operating condition, un- 
less operation under such conditions can be included in the stage characteristics. It 
should be noted that the simulation is not a valid representation of the compressor once 
the stability boundary is crossed. 


REMARKS AND CONCLUSIONS 

A high frequency analytical model for a jet engine compressor has been presented. 

The modeling technique was applied to the eight-stage axial-flow compressor of the 
J 85-13 engine. The stability boundary of the compressor simulation agrees closely with 
the stall line of the actual jet engine compressor. 

It was found, that it is necessary to include lumped volume representations for the 
gas dynamics of each compressor stage to achieve agreement between the simulation sta- 
bility boundary and the compressor stall line. 

Plots of the stage pressure characteristic for each compressor stage were presented 
for two compressor corrected speeds. The normal stage operating point and an operating 
point just prior to the stability boundary were indicated on each of the pressure character- 
istic curves. From these plots it was concluded that the stability boundary is a function 
of the combined nonlinear system gains and the volume dynamics. 

The compressor stability boundary is strongly influenced by the stage characteristics. 
The overall results are thus dependent on the accuracy of the stage characteristics, and 
the stall line prediction is dependent on the characteristics at the lower values of the flow 
parameter. The representation is one dimensional, and, thus, two-dimensional phenom- 
ena such as rotating stall and the effects of spatial distortion of the inlet flow field cannot 
be evaluated. If the technique proves successful in the prediction, the stall characteris- 
tics of other compressors, it will provide a useful tool for future compressor design pro- 
grams. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, July 10, 1969, 

720-07. 
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APPENDIX A 


SYMBOLS 

a ,,2 2 

A area, ft ; m 

a speed of sound, ft/ sec; m/sec 

C coefficient, rpm^/°R 2 ; rpm^/K^ 

C p specific heat at constant pressure, Btu/(lbm)(°R); J/(kg)(K) 

C y specific heat at constant volume, Btu/(lbm)(°R); J/(kg)(K) 

g gravitational constant, unit force/unit mass 
h enthalpy, Btu/lb m ; J/kg 

J mechanical equivalent of heat, ft-lbf/Btu; N-m/J 

K coefficient, (rpm)(see)/ft; (rpm)(sec)/m 

k bleed flow coefficient, °R/sec; K/sec 
L length, ft; m 

M Mach number 

N rotational speed, rpm 

P pressure, lb/ft^; N/m^ 

R universal gas constant, ft/°R; m/K 

r radius, in.; m 

T temperature, °R; K 

t time, sec 

U rotor speed, ft/sec; m/sec 

u internal energy, Btu/lb; J/kg 

V volume, ft^; m^ 

v velocity, ft/sec; m/sec 

v„ axial velocity, ft/sec; m/sec 

W weight flow, lb/sec; kg/sec 

a coefficient, 360 (in. )(sec)/(ft)(min); 30 (m)(sec)/(m)(min) 

/ 3 rotor air inlet angle, deg 
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y 

6 

0 

p 

$ 



Subscripts: 

b 

c 

n 

r 

s 


ratio of specific heats 

ratio of total pressure to standard atmospheric pressure nondimensional 

ratio of total temperature to standard atmospheric temperature 

weight density, lb/ft^; kg/m^ 

flow coefficient 

pressure coefficient 

temperature coefficient 

variable associated with stage bleed 

variable associated with the stage characteristic 

stage number designation 

reference state 

static condition 


t 


total condition 


v variable associated with the stage volume 

1, 2, 3, . , , , 8 stage numbers 

Superscript: 

- maximum value of quantity used for normalizing 
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APPENDIX B 


STEADY-STATE COMPRESSOR REPRESENTATION 


The analytical modeling of the steady-state compressor characteristics is based on 

the stage- stacking technique. This technique, which is discussed in reference 1, uses a 

pair of two-dimensional curves for each stage. The individual stage performance, de- 

P T 

scribed by a pressure coefficient and a temperature coefficient , is assumed to 
be a function of a flow coefficient <£. The coefficients are mathematically described by 


p 2gJ Ah. 
V = — — — E 


U 


T 2gJ Ah. 
i 

U 2 


$ = 


u 








(Bl) 


J.T- 

The ideal total enthalpy rise for the n stage can be related to the stage total pressure 
ratio from 


Ah tc, n C p T tv, n- 1 


p \ (r- 1 )/r 

tc, n \ _ 1 


L \V n- y 


(B2) 


The adiabatic enthalpy rise can be determined from 


Ah tc, n ~ C p( T tc,n " T tv, n-P 


and the stage wheel speed U from 


(B3) 


U = IS£ (B4) 

a 

The stage inlet conditions are associated with the conditions in the preceding volume and 
the variables are subscripted accordingly. The stage exit conditions are, however, as- 
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4* Vi 

sociated with the stage itself. Assuming that the n n compressor stage is preceded by 
the n - 1 volume, equations (Bl) can be combined with (B2) to (B4) to form the following 

VVi 

set of equations for the n ul stage 




P _ C n T tv, n-1 
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tc, n 
5 tv, n- 1/ 


,2/7 


^n 9 < T tc,.n " T tv, n-1^ 
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$ = -JL v 

n N zc, n 


(B5) 


where 


2gJC (a)' 

c„ = — — 2 — 


n 


2 2 
7r r 

n 


The flow coefficient $ is defined as the ratio of the axial flow velocity to a mean rotor 
velocity. Introducing the square root of the stage corrected temperature into the relation 
for $ given in equations (Bl), the flow coefficient becomes 


vJV® 
V 1W 


(B6) 


The axial velocity v„ can be computed from the flow, total pressure and total tempera- 

z 

ture at the stage inlet. The flow and the axial velocity are related by 


W 

Ps A 


Eliminating the density from (B7) with a state equation of the form 


P s Ps RT s 


(B7) 
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results in 


v 


z 


WR T s 
A P s 


(B8) 


or 


v 


z 


WR T t L | r - 1 

A P t V 2 



l/(r-i) 


Equation (B9) can be rewritten in terms of corrected parameters as 




l/(y-l) 


(B9) 


or 


WV0~ 

A6 


V© \ 2 


M 


2 \ - V (y ” i) 


p tr 


It can be shown from the energy equation that 

a t 2 >- R \ T t / 



thus, since 







(BIO) 


(Bll) 


(B12) 
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Combining equations (BIO) and (B12) results in 


A6 


\/© 



' v \ 2 1 

,V©/ 2gJC p T tr 




p tr 


(B13) 


It should be noted that the nonsubscripted velocity variable v introduced in equa- 
tion (Bll) is the magnitude of the velocity vector and not the axial velocity component 

v_. If the air flow approaches the rotor inlet at an angle 0 
z 


= V cos 0 
z 


and equation (B13) becomes 


WV0 _ v z 
A6 s/q" 




2 g JC p T tr cos ' P 


|i/(y-i) 


p tr 


For the n^ 1 


stage the flow parameter 


<3> n can be described by 



v zc, n 

N V®y,n-1 
V 0 v ,n-1 


(B14) 


(B15) 


(B16) 


and the corrected axial velocity v zc n /V® v n _i * s determined from solutions of equa- 
tion (B15) rewritten for n^ stage as 


W c, nV®v, n- 1 _ v zc, n 


^n 5 v , n-1 V®v, n- 1 


v 


zc, n 


v, n-l/ 2gJC p T tr cos 0 n J 


i/(y-i) 


P tr (B17) 
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APPENDIX C 


STAGE GAS DYNAMICS 


The gas dynamics of each compressor stage is assumed to occur in an equivalent 
stage volume. The stage volume is conceptually located just downstream of the stage 
pressure and temperature rise characteristics. The stage gas dynamics are represented 
by continuity, energy, and momentum balances across the volume. 

Assuming no work is performed in the volume and that the net heat input is negligible, 

1L. 

an energy balance can be written for the n 1 volume as 


cl 

^n ^sv, n u tv, n 


^ ^tc, n^c, n ” ^tv, n^e, n+1 


(Cl) 


Noting that u^ = C y T t and that h^ = C T^., equation (Cl) becomes 


^sv, n T tv, n^ “ y ( T tc, n^c, n ” T tv, n^c, n+P 


(C2) 


To determine the density of the gas in the stage volume, it is assumed that the rate of 
change of the static density is proportional to the net weight flow: 


^-(P B v,n>=7-(W C) n-W C!n+1 ) 


(C3) 


n 


The total pressure in the stage volume is calculated from an equation of state of the form 


2 \ 5 / 2 

P tv,n = R l 1 + 0 - 2I <nJ <P 


sv, n 


' tv, n' 


(C4) 


A momentum balance written across the stage and its associated volumes is used to 
provide an expression for the determination of weight flow. It is assumed that flow ac- 
celeration is caused by the difference between the static compressor head rise AP gc n 
and the head rise (P_ r „ - P t ) imposed by the adjacent volumes. Neglecting con- 
vective momentum terms, the above assumption yields 


_d 

dt 
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\_ s c , n 
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sv, n 


sv, n 


-!>] 


(C5) 
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The stage inlet conditions, however, are the conditions in the inlet volume and, thus, 
AP sc, n = P sc, n ' P sv, n- 1 and ec * uation < C5 ) becomes 


— (W n )=^(P 
dt c ’ n L. 


sc, n 


P ) 
sv, n' 


(C6) 


n 


Assuming the Mach number at the stage exit and in the volume are the same, equa- 
tion (C6) can be written as 


_d_ 

dt 


’ ( p tc,n 
Hi 


-7/2 


p tv,n> 


(C7) 


Equations (C2) to (C4) and (Cl) provide a representation for the gas dynamics of the 

iL. 

n m stage volume. In the computer simulation of these equations, the distinction between 
static and total pressures was dropped. This simplification resulted in neglecting the 
Mach number corrections in equations (C4) and (C7). The interstage Mach numbers en- 
countered in the compressor are less than 0. 5. The range of the corrections in equa- 
tions (C4) and (C7) are thus, 


/ ON 5 / 2 

1.0 < (l + 0.2 M'j <1.13 

/ o\ — 7/2 

0. 84 < (l + 0. 2 M^j <1.0 

The error introduced by this simplification is, thus, approximately equivalent to a 
15 percent change in the equation coefficients in the worst case. With the computer 
formulation used, the approximation is reflected in dynamic coefficients and does not 
alter the steady- state results. 

Figure 11 contains a block diagram of the computational scheme used for the dynam- 
ics of a typical stage. The coefficients used in the simulation are represented by num- 
bered circles and the neglected corrections by lettered dashed circles. Noting that the 
integrator inputs are nulled in the steady state, it can be shown that neglecting the Mach 
number corrections results in 
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T tv,n-1 



Figure 11. - Block diagram for stage computational scheme. 




in the steady state. Thus, the steady-state results are unaffected by the approximations 
in the volume dynamics. 
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